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Abstract.  The  Naval  Observatory  Flagstaff  Station  1-m  telescope  is 
evaluated  for  the  addition  of  adaptive  optics  capabilities  to  its  instrumen¬ 
tation  suite.  Zernike  decomposition  of  the  optical  system  based  on  phase 
diversity  measurements  shows  that  the  static  optical  aberrations  are 
small  enough  that  they  will  not  degrade  the  performance  of  the  deform¬ 
able  optical  element.  The  analysis  makes  use  of  annular  pupil  Zernike 
polynomial  reconstruction  of  the  wavefront  to  accommodate  the  large 
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1  Introduction 

The  Naval  Observatory  Flagstaff  Station  (NOFS)  1-m  tele¬ 
scope  was  the  last  telescope  to  be  designed  and  built  by  G. 
W.  Ritchey  and  is  a  Ritchey-Cretien  Cassegrain  optical  de¬ 
sign.  The  telescope  was  completed  in  1934  and  was  origi¬ 
nally  located  in  Washington,  D.C.  The  telescope  was 
moved  to  Flagstaff,  Arizona,  in  1955  to  take  advantage  of 
the  darker  skies.  This  telescope  has  been  the  mainstay  of 
work  at  the  Flagstaff  Station  and,  in  keeping  with  this  role, 
it  is  now  being  evaluated  for  upgrading  to  high  angular 
resolution  work  by  adding  an  adaptive  optics  system. 

Other  telescopes  built  by  Ritchey  include  the  60-  and 
100-in.  telescopes  at  the  Mount  Wilson  Observatory,  which 
have  both  had  adaptive  optics  systems  retrofitted  to  them  at 
various  times  in  the  past.  These  telescopes  are  conventional 
Cassegrain  telescopes  with  obscurations  less  than  about 
28%.  The  unique  feature  of  the  NOFS  1-m  telescope  that 
makes  the  analysis  more  interesting  is  that  the  secondary 
mirror  is  particularly  large,  having  an  approximately  45% 
obscuration,  which  was  designed  to  provide  an  unvignetted 
field  of  view  of  greater  than  a  degree.  While  it  has  been 
known  for  many  years  that  the  use  of  filled  circular  Zernike 
polynomials  on  telescopes  with  large  obscurations  can  lead 
to  significant  errors,  a  recent  literature  search  has  shown 
that  few  papers  exist  on  the  use  of  the  more  appropriate 
annular  Zernike  polynomials1  for  analysis. 


This  paper  presents  an  analysis  of  the  existing  optical 
system  of  the  NOFS  1-m  telescope  to  evaluate  whether  the 
static  aberrations  in  the  existing  telescope  optics  would  re¬ 
duce  the  effectiveness  of  an  adaptive  optics  system.  These 
measurements  are  based  on  wavefront  measurements  using 
the  technique  of  phase  diversity2  with  Zernike  polynomial 
decomposition  used  to  evaluate  the  telescope’s  wavefront 
quality.  An  annular  Zernike  polynomials  expansion  was 
chosen  for  this  analysis  because  of  the  telescope’s  large 
obscuration  and  the  results  are  compared  with  using  con¬ 
ventional  filled  circular  Zernike  polynomials. 

2  Telescope  Wavefront  Analysis 

The  “as  built”  prescription  for  the  NOFS  1-m  telescope 
was  evaluated  in  Zemax  and  was  used  to  simulate  the 
wavefront  characteristics  of  the  telescope.  This  analysis 
showed  that  the  peak-to-valley  wavefront  error  was  better 
than  \/5.  In  these  calculations,  the  distance  from  the  sec¬ 
ondary  mirror  to  the  focal  plane  was  optimized  for  the  best 
focused  spot  size  for  this  simulation.  It  was  also  clear  from 
the  Zemax®  model  that  the  telescope’s  large  obscuration 
represented  a  significant  deviation  from  the  assumption  of 
“no  obscuration”  for  the  use  of  filled  circular  Zernike  poly¬ 
nomials.  To  analyze  the  wavefront,  the  technique  of  phase 
diversity  was  selected  and  the  results  from  filled  circular 
and  annular  pupil  Zernike  polynomial  decomposition  are 
compared  in  the  evaluation  of  this  wavefront. 
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Fig.  1  Raw  out-of-focus  images  taken  at  the  NOFS  1-m  telescope. 


Stellar  images  were  obtained  from  the  NOFS  1-m  tele¬ 
scope,  recorded  both  inside  and  outside  focus,  as  shown  in 
Fig.  1,  for  evaluation  using  the  phase  diversity  technique.2 
The  objective  was  to  analyze  the  wavefront  that  is  deliv¬ 
ered  from  the  telescope  to  determine  whether  or  not  there 
were  significant  optical  aberrations  present  in  the  telescope. 
For  this  analysis,  we  considered  the  transport  equation  for 
intensity  and  phase,3'4  as  shown  in  Eq.  (1). 

dl 

-k  —  =  V-(/V</>)  =  V/V<£+/V2<£,  (1) 

dz 

where  I  is  the  intensity,  cp  is  the  phase,  k  =  2rr/\,  and  z  is 
the  propagation  axis.  The  derivative  along  the  optical  axis 
is  approximated  with  a  finite  difference.  In  the  most  com¬ 
mon  case,  the  two  intensity  maps  recorded  from  both  sides 
of  the  Gaussian  focus  are  used  to  generate  All  Az .  This  is 
used  to  solve  Eq.  (1).  Assuming  a  nearly  constant  intensity 
map  in  the  plane  of  observation  reduces  Eq.  (1)  to 

dl 

-k-=/V2</>,  (2) 

dz 

where  Eq.  (2)  is  the  usual  form  of  the  curvature  sensing 
equation.-4 

The  images  in  Fig.  1  were  obtained  at  the  telescope  by 
recording  the  out-of-focus  images  from  a  bright  star.  These 
images  were  taken  at  comparable  distances  from  the  best 
focus  position.  The  images  were  recorded  using  an  expo¬ 
sure  time  considerably  longer  than  the  “coherence  time  of 
the  atmosphere”  to  average  out  the  contributions  of  the 
atmospheric  turbulence.  Typical  values  for  the  “coherence 
time  of  the  atmosphere”  at  a  good  site  range  around  10  m. 

Figure  2  shows  the  difference  of  the  two  intensity  im¬ 
ages,  corresponding  to  AH  A z  on  the  left  and  the  phase 
image  calculated  from  it  using  Eq.  (2)  on  the  right.  In  Fig. 
2  (right)  only  two  spiders  can  be  seen  after  the  conversion 
from  an  intensity  image  as  compared  to  the  four  that  are 
easily  seen  in  the  raw  images  in  Fig.  1.  This  can  be  attrib¬ 
uted  to  the  presence  of  astigmatism  enhancing  the  second¬ 
ary  spider  support’s  footprint  along  one  diagonal  and  reduc¬ 
ing  the  phase  contrast  of  the  opposite  diagonal. 

It  is  convenient  to  decompose  the  phase  wavefront  into 
an  orthonormal  basis  set  (Zernike  modes)  to  evaluate  the 
wavefront.  For  sake  of  simplicity  in  plotting  the  decompo¬ 
sition  of  the  wavefronts  we  count  the  Zemike  modes  in  a 
sequential  manner  instead  of  using  the  customary  2-D  in¬ 


Fig.  2  Left:  difference  of  the  raw  out-of-focus  images  of  Fig.  1,  an 
approximation  of  the  dlldz  component  of  Eqs.  (1)  and  (2);  Right: 
phase  map  calculated  from  the  left  image  using  Eq.  (2).  Notice  that 
there  is  an  enhancement  of  the  spiders  in  one  direction  while  they 
are  lost  in  the  orthogonal  direction. 


dexing;  the  indexing  order  is  shown  in  Fig.  3,  beginning  in 
the  lower  left  corner  of  the  plot,  where  n  and  \m\  are  the 
radial  and  azimuthal  indices,  respectively. 

Traditionally,  filled  circular  Zernike  polynomials  are 
used  for  this  type  of  decomposition  and  analysis.  In  the 
case  of  the  NOFS  1-m  telescope,  the  large  obscuration 
would  suggest  that  annular  Zernike  polynomials  would  be 
more  appropriate,  since  these  polynomials  are  an  orthonor¬ 
mal  basis  set  for  an  annulus.  To  examine  the  difference 
between  these  two  variants  of  Zernike  analysis  the  mea¬ 
sured  wavefronts  are  decomposed  into  the  amplitudes  of 
the  coefficients  of  the  modes,  which  are  plotted  as  bar 
charts  in  Figs.  4  and  5.  These  values  were  calculated  using 
custom  software  built  in  IDL  (Ref.  3)  and  the  mode  num¬ 
bering  used  is  that  shown  in  Fig.  3. 

The  measured  wavefronts  are  composed  of  optical  aber¬ 
rations  from  the  telescope  misalignments  and  static  aberra¬ 
tions  in  the  optics  themselves,  while  atmospherically  in¬ 
duced  aberrations  are  minimized  by  the  averaging  effects  of 
the  exposure  time.  The  aberrations  induced  by  alignment 
errors  cannot  be  easily  separated  from  the  static  optical 
aberrations  in  this  case  and  are  assumed  to  be  much  smaller 
than  the  atmospherically  induced  ones.  Comparing  these 
two  figures  it  can  be  seen  that  the  filled  circular  case  shows 


M 


Fig.  3  Single  index  ordering  of  the  2-D  Zernike  array  used  com¬ 
pared  to  the  customary  2-D  indexing.  The  sequence  begins  in  the 
lower  left  of  the  plot. 
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Fig.  4  Magnitude  plot  of  the  coefficients  for  the  filled  circular  Zernike  polynomial  modes. 


Zernike  expansion  (annulor  pupil} 


Fig.  5  Magnitude  plot  of  the  coefficients  for  the  annular  pupil  Zernike  polynomial  modes. 
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Table  1  Comparison  of  filled  circular  versus  annular  pupil  Zernike  polynomials. 


Spherical  aberration 
Coma 

Astigmatism 

Defocus 


A4(6P4-6P2+1) 

A3(3P3-2P)cos  0 

A|p2  cos  20 
A2(2P2-1) 


A4' 

V5(1  -e2)2 


[6p4— 6(1  +E2)p2+1  +4e2+e4] 


2  1  +e2+e4  \ 

3  TTe2  Pjcos  ® 


a2 i  (2  cos  20-1  )p2 

Az  ( 1  +e2+e4)1/2 

^o’[p2-(1+e2)/2] 


several  different  effects.  The  relative  importance  of  the 
modes  between  the  two  Zernike  expansions  is  similar,  that 
is,  the  fractional  weight  of  the  amplitude  of  the  modes  are 
comparable,  however,  it  is  significant  that  in  many  of  the 
modes,  the  associated  sign  has  reversed.  As  examples,  in 
the  first  five  terms,  modes  1  and  4  show  opposite  signs 
when  the  two  plots  are  compared.  The  effect  of  using  filled 
circular  versus  annular  pupil  Zernike  modes  is  discussed  in 
the  next  section. 

3  Discussion 

The  large  obscuration  of  the  NOFS  1-m  telescope  was  the 
original  motivation  for  examining  the  use  of  annular  pupil 
Zernike  polynomials1  in  the  analysis.  Table  1  compares  the 
form  of  selected  filled  circular  Zernike  polynomials  with 
the  annular  pupil  Zernike  polynomials.  In  the  table,  Aj!m|  is 
the  coefficient  of  the  mode,  p  is  the  radial  component,  0  is 
the  angular  component,  and  e  is  the  fractional  diameter  of 
the  obscuration.  In  Table  1  we  identify  the  Zernike  mode 
with  the  correspondent  balanced  aberration  as  is  customar¬ 
ily  done  (see,  for  example.  Ref.  5). 

The  effect  of  the  use  of  the  annular  pupil  Zernike  repre¬ 
sentation  is  compared  with  the  filled  circular  Zernike  for 
the  two  cases  of  spherical  aberration  and  astigmatism 
shown  in  Fig.  6.  In  these  two  cases,  the  change  in  the 
magnitude  of  the  polynomial  can  be  clearly  seen,  however, 
the  astigmatism  plot  shows  how  the  sign  change  for  the 
mode  occurs.  In  adaptive  optics  systems,  this  is  not  a  par¬ 


ticularly  important  effect  as  typically  the  wavefront  analy¬ 
sis  does  not  use  the  Zernike  modes  explicitly  but  rather  sets 
the  deformable  element  to  null  out  the  wavefront  error. 
However,  this  would  be  important  in  the  case  of  determin¬ 
ing  the  prescription  of  the  optic. 

Figure  7  shows  the  wavefront  that  is  formed  by  recom¬ 
bining  40  annular  pupil  Zernike  modes  from  the  phase  di¬ 
versity  data  following  Eq.  (2).  In  this  image,  several  low- 
order  features  can  be  seen  including  evidence  of  a  mixture 
of  trifold  coma  and  astigmatism.  This  aberration  is  also 
evident  in  the  amplitudes  of  these  modes  shown  in  Fig.  5. 
The  root  mean  square  (rms)  value  of  the  difference  between 
the  phase  map  shown  in  Fig.  2  (right)  and  the  two  Zernike 
expansions,  full  pupil  and  annular  pupil  is  rmsann=0.05  and 
rmSfuu=0.08,  respectively,  both  with  a  standard  deviation 
of  0.002.  In  either  of  the  Zernike  models,  the  rms  wave- 
front  errors  are  quite  small,  which  enables  their  correction 
using  deformable  mirrors  or  spatial  light  modulators  to  cor¬ 
rect  these  terms  easily.  However,  should  there  be  a  desire  to 
retouch  the  optical  components  themselves,  how  the  coef¬ 
ficients  are  distributed  across  the  two  reflective  elements  of 
the  telescope  should  be  determined. 

Based  on  this  analysis,  the  amount  of  static  aberration5 
present  in  this  telescope  is  roughly  A./4  in  magnitude  and 
mostly  concentrated  in  the  low-order  aberrations,  as  seen  in 
Fig.  5.  For  most  adaptive  optics  systems,  static  errors  of 
this  order  can  be  easily  corrected  without  compromising  the 
dynamic  range  of  the  system.  Much  of  the  concern  for  this 


Fig.  6  Cross-cuts  of  the  Zernike  mode  for  astigmatism  (left)  and  spherical  aberration  (right)  for  filled 
pupil  (solid  line)  and  annular  pupil  (dashed  line). 
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Fig.  7  Reconstruction  of  the  wavefront  using  a  40-term  expansion 
of  the  annular  pupil  Zernike  fit  of  the  data. 

telescope  was  due  to  the  large  secondary  obscuration  and 
the  effect  that  this  may  have  on  the  number  of  illuminated 
actuators  in  the  pupil  available  for  correction.  Since  many 
of  the  active  elements  of  the  wavefront  corrector  would  be 
hidden  in  the  shadow  of  the  secondary  and  thus  not  illumi¬ 
nated  and  available  for  use,  the  effectiveness  of  the  wave- 
front  corrector  could  be  significantly  reduced. 

4  Summary 

A  phase  diversity  analysis  of  the  NOFS  1-m  telescope  was 
carried  out.  The  results  show  that  the  rms  wavefront  error, 
which  we  have  assumed  as  a  static  error  in  the  optics,  is 
small  and  easily  correctable  using  a  wide  variety  of  stan¬ 
dard  wavefront-correcting  elements.  Using  such  optical  el¬ 
ements  to  correct  the  static  aberrations  will  leave  sufficient 
residual  throw  in  the  deformable  element  to  compensate  for 
atmospheric  seeing. 
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